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The neutrino propagation in neutron stars is studied in the framework of the linear response 
method. The medium effects are treated in the non-relativistic Brueckner-Hartree-Fock approach 
either in the mean- field approximation or in the RPA. The residual interaction is expressed in terms 
of the Landau parameters extracted from the equation of state of spin- and isospin-polarized nuclear 
matter. The Brueckner theory including three-body forces is used for determining the equation of 
state. Numerical predictions for the response function of nuclear matter in /3-equilibrium and the 
neutrino mean free path are presented in a range of baryonic densities and temperatures. The main 
results are a dominance of the charge-exchange component over the scattering component and an 
enhancement of the neutrino mean free path induced by nuclear correlations. 

PACS numbers: 26.60.+C, 26.50.+X 
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I. INTRODUCTION 

The interaction of neutrinos with baryons has been 
mostly studied in connection with the stability of nuclei, 
but it also plays a crucial role in the thermal evolution 
of supernovae and protoneutron stars, where the nuclear 
medium exhibits quite distinctive physical features re- 
lated to density, temperature and chemical composition. 
A large effort has been devoted to study the produc- 
tion of neutrinos via direct or modified URCA processes 
1] and, more recently, via the bremsstrahlung of nucle- 
ons in the strong magnetic field of neutron starsQ- In 
the past simulation codes have incorporated the propaga- 
tion of neutrinos in neutron matter mostly by using the 
Fermi-gas model, whereas the effects of correlations did 
not receive the due consideration except in a few cases 
El 01 ■ In the recent years calculations of correlation 
effects have been performed, including phenomenologi- 
cal approaches with Skyrme forces |fj, [D, |8|, microscopic 
Brueckner theory ^ IM P |Io | and relativistic mean field 
theory [ill 13 ES GjJjllllM El On the other hand, 
in order to reliably describe the correlations of nuclear 
matter in extreme conditions, one needs a well devel- 
oped microscopic many-body theory, since not enough 
experimental constraints exist so far in a so wide range 
of density, temperature and isospin asymmetry which are 
supposed to occur in the present neutron-star models. 

In this paper we present a study of the nuclear re- 
sponse function to weak interaction and the effects of 
short- and long-range correlations on the neutrino trans- 
port in neutron stars. It is based on the equation of 
state (EOS) predicted by the Brueckner theory approach 
including relativistic effects and nucleonic resonances. 
Recently this approach has made a step forward in re- 
producing the empirical saturation properties of nuclear 
matter 19, 20]. The particle-hole NN residual interaction 
is extracted from the Brueckner theory and cast in terms 
of the Landau parameters. The response function to the 



neutrino propagation is calculated first in the Brueckner- 
Hartree-Fock (BHF) limit, and then in the random phase 
approximation (RPA) in neutron matter in /3-equilibrium 
with protons and electrons for various conditions of neu- 
tron density and temperature. The neutrino mean free 
path (MFP) is also calculated taking into account the 
(y e , v' e ) scattering from nucleons and the (y e , e~) charge- 
exchange process. In this work we discuss only the v e 
case which is the most important one for the MFP issue. 



II. NUCLEAR MATTER IN THE BHF 
APPROXIMATION 

A neutron star, in the most simplified model, is made 
of neutrons, protons and electrons, and the relative abun- 
dance is controlled by the /3-equilibrium condition (the 
effect of muons is negligible). Under the condition of 
charge neutrality, the proton fraction Y p — Z/A is driven 
by the value of the symmetry energy at a given total 
baryonic density. Here, we will always consider nuclear 
matter in /3-cquilibrium unless explicitly stated. 

Our framework is the Brucckner-Bethe-Goldstonc 
(BBG) approach, where the perturbative expansion of 
the total energy per particle Ea can be cast according 
to the number of hole lines. The lowest order defines the 
Brueckner-Hartree-Fock approximation, which exhibits a 
satisfactory convergence provided the continuous choice 
for the mean field is adopted ■ 

In the recent years a remarkable step forward in the 
reproduction of the empirical saturation properties has 
been performed including three-body forces (3BF) in the 
BBG theory [Hill. Their effect in fact is twofold: on 
one hand, it brings the saturation density of nuclear mat- 
ter close to the empirical value; on the other hand, it pro- 
vides the high-density strong repulsive components. As 
a consequence the BHF approximation has reached the 
level of a consistent description of a strongly correlated 
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Fermi system in the region of saturation density as well 
as in high density regions. 

A microscopic derivation of the full set of I = Landau 
parameters can be done, based on the BHF approxima- 
tion, from the calculation of the energy per particle of 
nuclear matter in different states of spin and isospin po- 
larization. Starting from unpolarized symmetric nuclear 
matter at a given density p, we may extend the calcula- 
tion of Ea to isospin density variations Sp T = (p T — p)/p 
and spin-isospin density variations 8a T = (pi — p\.)/p T . 
The response of the system is related to the second 
derivatives of the expansion of Ea, 

E A (p,5p n ,5p p ) = E A (p) + - ^2$ t ,t>SPtSPt>, (1) 

tt' 

E A {p, 5<j n , 5<t p ) = Ea{p) + - r T . Tl 5<j T 5<r T >. (2) 



These derivatives give the / 
the residual NN interaction, 
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where iV(0) is the level density at the Fermi surface. No- 
tice that our definitions Eqs.(3)-(6) differ by a factor of 2 
from that of other authors, for instance ■ The Landau 
parameters are constrained by robust physical observ- 
ables. Fq is in fact constrained by the compression mod- 
ulus, whose value should range in the interval 230-250 
MeV according to the experimental determination of the 
monopole giant resonance energy centroid. Fq is related 
to the symmetry energy, which at the saturation density 
is reported about 30 MeV in the Bethe-Weiszacker mass 
formula. Go is related to the spin modes, which actually 
are hardly observed in nuclei. So far experimental infor- 
mation on Go is not enough since spin resonances have 
only been observed with too small strength compared to 
other collective modes [22J. Finally G' is constrained by 
the Gamow- Teller giant resonance. A value of 1.2 at the 
saturation point has been determined with high precision 
from the experimental excitation energy of the Gamow- 
Teller resonance on 90 Ni |2_3]. The BHF prediction of the 
Landau parameters is reported in Fig. 1 . They have been 
obtained from a BHF calculation with the Argonne AVis 
[2^ as two-body force (2BF) and a microscopic three- 
body force which corrects the BHF approximation by 
relativistic effects and nucleonic virtual excitationsp5|. 

The prediction of the Landau parameters for densities 
other than the nuclear density is of great interest in the 
study of neutron stars. In connection with the strong 
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FIG. 1: Landau parameters of symmetric nuclear matter with 
2BF AVis (right) and the same plus 3BF (left) 



magnetic helds observed in neutron stars some authors 
have studied the magnetic susceptibility in neutron mat- 
ter and found that Go reduces the susceptibility of the 
degenerate neutron gas|2^, |27l IH, |2^| . This reduction 
is amplified at high de nsity when including 3BF either in 
Brueckner calculations [23 or in Monte Carlo many-body 
simulations pfj . 

In this work it is assumed that the Landau parame- 
ters do not change appreciably with temperature, which 
is certainly a good approximation below the liquid-to- 
gas phase transition (T c w 18MeV). This point will be 
discussed in more details in section IV(A). 



III. INTERACTION OF NEUTRINOS WITH 
MATTER 

During their propagation in neutron-star matter neu- 
trinos experience collisions with nucleons via weak cou- 
pling with the nucleon neutral currents j M = t\> ' T 7 '^(c£ — 
Ca"/c,)ipT- Neutrinos can also disappear in the charge- 
exchange process v e +n — > e~ +p and we also consider the 
coupling with charged currents = ip^^gv - SVCs)^ 
giving rise to neutrino absorption. The vector and ax- 
ial coupling constants for neutral currents are c£ and c„, 
and for charged currents g v and g a (a complete list of the 
constants is given in Ref.|3)- 

The neutrino MFP A is derived from the transport 
equation 

£ = W fi (k - fc')(l - n % ,) + W fi (k' - k)n %n (7) 



where is the occupation number of neutrinos and Wfi 
is the transition probability corresponding to the individ- 
ual processes under consideration. 

In the non-relativistic limit the transition rate can be 
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written[! HQ 



Wfi = G 2 W [(1 + cos 6) W v + (3 - cos 0)W a ] 



(8) 



where G w is the weak coupling constant, 9 is the scat- 
tering angle and W v (W a ) is the vector (axial) term. 
The transition rates are expressed in terms of the nu- 
clear structure functions S^ T a . In the case of nucleon 
scattering the relation is 

W v<a = ^ ' c v,a c v,a^v,a J (9) 
t,t' 

< - - l ~Yz^ lm ^' (10) 

For the charge- exchange process the relation is 

W»,a = ffL^"a> (11) 

e TT ' = -1 - TmY"' a ("12") 

tt 1 _ e (qo+S^)/kT LLLL ^P,n- \ L *> 

The quantity q is the energy transfer and Sfi is the 
shift between the neutron and proton chemical poten- 
tials. From the previous equations we see that the struc- 
ture functions are, in turn, related to the response func- 
tions x which will be discussed in the next section. 



IV. NUCLEAR RESPONSE FUNCTIONS 

The major part of the work needed to obtain the neu- 
trino MFP is in the calculation of the structure functions. 
These calculations can be performed at different levels 
of approximation beyond the free Fermi gas model. We 
will consider first the BHF approximation, which already 
takes into account the strong short-range correlations via 
the mass renormalization and the depletion of the Fermi 
surface; second, the particle-hole interaction will be taken 
into account in the RPA approximation. Since the ring 
diagram summation is almost a prohibitive task with the 
Brueckner G-matrix, we may approximate, in the low- 
frequency limit bjKkvp) <C 1, the G-matrix residual in- 
teraction by a Landau-Migdal form whose parameters are 
extracted from the BHF theory as described in the pre- 
vious section. 



fulfilled. Once the approximation for £ has been settled, 
the single-particle propagator takes the form 

G^fow) = w -|--£ T (p, w )+ e J, « ^f0"-#, (13) 



where t r v is the single-particle energy 
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and €p is the Fermi energy. The factor = (1 — 

^r)^ie p i s l ne quasiparticle strength associated to the 
depletion of the occupation probability of the single- 
particle level e p . In the present study the rearrangement 
term will be neglected, which is not a too severe approxi- 
mation since it does not affect the EOS at zero tempera- 
ture and, consequently, the calculation of the Landau pa- 
rameters . At finite temperature it cannot be absolutely 
neglected, since the single-particle spectrum determines 
the thermal Fermi functions and hence the EOS itself. 

The neutrino propagation is studied in the frame- 
work of the linear response theory. Let us first discuss 
the BHF limit where there are no correlations coming 
from the residual interaction. The dynamical structure 
functions are related to the imaginary part of the re- 
tarded polarization or response function (<f, qo) = 

(Rc + isign(go)Im) n^,(g,go), iK ) being the polariza- 
tion function 



<},(&«to) = -2t J 
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where the propagators are calculated according to 
Eqs. (13-14). One easily finds 
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According to the previous discussion the thermal occu- 
pation numbers 
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A. BHF response function 



are calculated from the single-particle spectrum frozen at 
T = 0. 



In the BBG approach the expansion of the self-energy 
E T (p, w)can also be cast as a hole line expansion, but 
the lowest order in the self-energy expansion does not 
correspond to the lowest order in the expansion of the 
energy per particle. According to the Landau definition 
of the quasiparticle energy, the BHF approximation for 
the self-energy should include the expansion up to the 
second order, the latter being the so called rearrangement 
term j^J. Doing so, the Hugenholtz-Van Hove theorem is 



B. RPA response function 

In this subsection we study the effect of the resid- 
ual particle-hole interaction. As already anticipated, the 
residual interaction will be described in terms of the I = 
Landau parameters discussed earlier. This makes it eas- 
ier to calculate the response function X ^ (<T; Qo) m the 
RPA limit. 
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Let us consider first the nucleon scattering. In this case 
the RPA response function is given by the Bethe-Salpeter 
equation in the 2x2 isospin space 



Xs = X 



(0) 



X 



(0) L (S) 



Xs, 



(18) 



where only the diagonal elements x-ril come into play and 
2/1°*./ = F2^i and lS;L = G" , are the matrix elements 
of the particle-hole residual interaction expressed by the 
Landau parameters (see Eqs.(3)-(6)). The S = (S = 
1) component is the vector (axial) part of the response 
function. 

In the case of charge-exchange processes the Bethe- 
Salpeter equation is a one-dimensional equation 



Xs 



y (o) 
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(19) 



where = Fq, AfW = G' , and X pn is the off-diagonal 
matrix element of the BHF response function. 

The structure functions are essentially given by the 
imaginary part of the response functions, according to 
Eqs.(10) and (12). 



V. RESULTS AND DISCUSSION 

The input data, namely the BHF mean field for the 
response functions and the Landau parameters for 
the RPA response functions, are obtained from a non- 
relativistic BHF calculation of nuclear matter with the 
two-body ^4Vi 8 [24| and the microscopic 3BF 25]. The 
temperature is set to zero, which is actually a good ap- 
proximation only below the critical point, as already 
mentioned. The proton fraction is calculated from the 
symmetry energy extracted in the calculation of isospin- 
asymmetric nuclear matter. For p = 0.34 fm~ 3 , Y p — 
0.077 without the 3BF and Y p = 0.167 with the 3BF. 
Thus, the 3BF has a sizable effect on Y p . 

The unperturbed (BHF) response function, as de- 
scribed by Eq.(16), is plotted in Fig. 2 for different tem- 
peratures and different processes. As it can be seen in 
the left panels (imaginary parts), the Ix^^nn^ 1 )! is al- 
ways greater than |x (PP _1 )|- This effect, which arises 
from the small proton fraction, i.e., large asymmetry pa- 
rameter 
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in neutron stars and smaller Z£ compared to from 
BHF calculations, gives larger scattering probability for 
v + n — ► v + n and, consequently, shorter neutrino MFP. 
The response function related to scattering contributes 
to the neutrino MFP essentially for qo > whereas the 
one related to absorption contributes for qo > (p p — fin). 
It also shows that the response function only depends 
softly on the temperature, at least in the T = Brueck- 
ner approximation, since the temperature appears only 
through n r k . 
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FIG. 2: The uncorrelated response functions for (v, v') scat- 
tering on neutrons (nn" 1 , solid), on protons (pp" 1 , dash) 
and neutrino absorption (pn _1 , dot). The imaginary part (left 
panel) and the real part (right panel) are plotted as functions 
of energy transfer. The p, /3 and momentum transfer q are 
fixed at 0.34 fm -3 , 0.668, 1.0 fm _1 , respectively. The upper 
and lower panels correspond to T=0 and T=20 MeV. 



Fig. 3 shows the imaginary part of the RPA response 
function as calculated from the Landau parameters with 
3BF included. Here, we show the vector (left side) and 
axial (right side) components for the three isospin chan- 
nels in consideration. The comparison emphasizes the 
dominance in both 5 = and S = 1 channels of the ab- 
sorption component (see insets). This effect is due to the 
q threshold (compare Eq.(10) and Eq.(12)) associated 
with the charge-exchange process as already discussed 
above. 

The resonant structure of all curves is due to the col- 
lective excitations of the medium: zero sound in the case 
of the vector component and spin waves in protons and 
neutrons in the case of the axial component. 

With the above response functions the neutrino MFP 
can be evaluated. The calculated MFP are shown as 
functions of the density in Fig. 4 for neutral current pro- 
cess on the left side and charged current process on the 
right side. T and E v are fixed to 20 MeV and 40 MeV, 
respectively, and (3 is fixed according to /3-equilibrium 
for each density and different approximations. In or- 
der to see the different effects from the scattering pro- 
cess and absorption process, we divide the neutrino MFP 
into two parts: A scat t. for scattering process (left panel) 
and A a bsorp. for neutrino absorption process (right panel). 
Then, the final neutrino MFP is A = (A"^ + A^ )~\ 
Comparing the left and right panels, we get the first mes- 
sage that the neutrino MFP caused by the neutrino ab- 
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FIG. 3: Imaginary parts of RPA response functions. The vec- 
tor components (left panel) and the axial-vector components 
(right panel) are plotted as a function of energy transfer for 
fixed nucleonic density, temperature and momentum transfer. 
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FIG. 4: Neutrino MFP for scattering (left panel) and absorp- 
tion (right panel). Solid (dashed) lines: RPA with (without) 
3BF; dotted (dash-dotted) lines: BHF with (without) 3BF; 
dash-dot-dotted lines: free Fermi-gas model. The tempera- 
ture T and neutrino energy E v are fixed to 20 MeV and 40 
MeV, respectively. 



sorption is much smaller than that caused by neutrino 
scattering. This is a well-known effect which is due to 
two reasons, as noticed by Reddy et al.@. First, the 
absorption transition probability is four times the scat- 
tering transition probability; second, the absorption rate 
is controlled by the electron chemical potential while the 
scattering is controlled by the neutrino chemical poten- 
tial. The inclusion of RPA correlations in the structure 
functions has the effect of reducing the phase space avail- 
able for scattering and the MFP turns out to be substan- 
tially increased compared to the free Fermi-gas model 
and also to the BHF model, the latter only at high den- 
sity. According to Ref . |2(J , the proton fraction in neutron 
stars increases dramatically at higher densities when the 
3BF is included. This affects the behavior of the neu- 
trino MFP in both scattering and absorption processes: 
the neutrino MFP increases with increasing density, espe- 
cially for the neutrino absorption process. This is related 
to the fact that with decreasing (3, the proton Fermi en- 
ergy increases and then the reaction v e + n — * e~ +p has 
a smaller cross section. 

Since the temperature in the medium drops out rapidly 
in the protoneutron star, it is interesting to show the re- 
lation between neutrino MFP and temperature. In Fig. 5 
(left) we plot the neutrino MFP as a function of temper- 
ature with p, (3 and E v fixed at 0.34 fm~ 3 , 0.668 and 40 
MeV, respectively. One can see the different behaviour 
of the scattering and absorption MFP when T increases. 
For the scattering process, A decreases by a factor 30 
from T=5 MeV to T = 60 MeV while for the neutrino 
absorption process this factor is only 2.5 . In the right 
side of Fig. 5 the E v dependence is also depicted. We just 
notice that the ratio between scattering and absorption 
MFP is almost constant. 
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FIG. 5: Neutrino MFP vs. temperature (left) and vs. neu- 
trino energy (right). The solid line is for neutrino absorption 
and dashed line for neutrino scattering. The values of p, (3, 
are 0.34 fm" 3 , 0.668, respectively. 



VI. COMPARISON WITH OTHER 
APPROACHES AND CONCLUSIONS 

The nuclear response function and the neutrino MFP 
in (3- equilibrium nuclear matter have been calculated 
on the one hand in the BHF mean field approximation, 
and on the other hand in the RPA using a Landau- 
Migdal residual interaction to evaluate the effects of long- 
range correlations. Mean fields and residual interaction 
have been obtained within the non-relativistic Brueckner- 
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Hartree-Fock theory. The two-body NN potential has 
been supplemented by a microscopic three-body poten- 
tial, which takes into account relativistic effects and also 
nucleonic excitations. From the present calculations one 
may conclude that the reduction of the neutrino MFP 
in the limit of the BHF mean field is much less pro- 
nounced than in the Fermi-gas model when going from 
low to high density .This reveals the important role of 
the short-range correlations incorporated in the BHF nu- 
cleon propagators. In addition, the effects of the residual 
particle-hole interaction have also a deep influence on 
the neutrino transport in /3-equilibrium nuclear matter, 
resulting in a sizeable enhancement of the neutrino MFP. 
This result confirms the conclusion of a previous calcula- 
tion of neutrino MFP in pure neutron matter @ . At high 
density this effect is magnified in the absorption channel 
by the intense repulsive three-body force. 

Since the interaction is poorly known at densities of 
interest in neutron star physics, the comparison with 
other approaches is difficult and leads often to contradic- 
tory predictions. Among the non-relativistic calculations 
we should quote the pioneering results by Burrows and 
Sawyer In spite of using a different set of Landau 

parameters (in particular Go = 0) the main features of 
neutrino MFP they predict are in good agreement with 
ours. In particular, they find an increase of A„ by up 
to one order of magnitude with increasing density as an 
overall effect of medium correlations, and a similar tem- 
perature dependence. In the recent work of Margueron 
et al.|l(]| the emphasis is put on temperature effects in 
a calculation of MFP based on BHF and RPA but the 
3BF is not treated and the important contribution to 



the MFP due to neutrino absorption is not calculated. 
It is also found that the scattering MFP including RPA 
correlations does increase with density, similarly to our 
present results. 

Non relativistic phenomenological approaches 0, 
manifest an opposite trend: the RPA neutrino MFP is 
quenched by the density increase. In fact, the residual in- 
teraction in the spin channel instead of slowly increasing 
at high density, decreases and finally becomes attractive. 
This explains also the prediction of a ferromagnetic tran- 
sition in phenomenological approaches. 

Relativistic approaches lead to an enhancement of the 
neutrino MFP when including the interaction effects in 
RPA even if some disagreement still exists among differ- 
ent calculations. In one case the reduction of the 
scattering rate is estimated about 10% to 15%, in others 
[rllT^ the effect is much stronger. 

A general conclusion is that the neutrino MFP is very 
sensitive to the EOS of nuclear matter in all aspects and 
it must be calculated as much as possible consistently 
with the EOS. In this respect we should mention a serious 
drawback of most calculations, which is the temperature 
dependence of the effective interaction. A typical value 
of temperature of about 20 MeV maybe is large for ne- 
glecting the effects of thermal excitations of the medium 
on the in-medium interaction, for instance the G-matrix 
in the Brueckner approach. 

In conclusion, our calculation predicts a quenching of 
the opacity of nuclear matter to the neutrino propagation 
and then provides the quantitative basis to estimate the 
energy release in the supernovae explosions and the rapid 
cooling of a protoneutron star. 
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